Evidence for the Strong Dominance of Proton-Neutron Correlations in Nuclei 
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Abstract: We analyze recent data from high-momentum-transfer (p,pp) and (p,ppn) reactions 
on Carbon. For this analysis, the two-nucleon short-range correlation (NN-SRC) model for backward 
nucleon emission is extended to include the motion of the NN-pair in the mean field. The model is 
found to describe major characteristics of the data. Our analysis demonstrates that the removal of a 
proton from the nucleus with initial momentum 275 — 550 MeV/c is 92l* g % of the time accompanied 
by the emission of a correlated neutron that carries momentum roughly equal and opposite to the 
initial proton momentum. Within the NN-SRC dominance assumption the data indicate that the 
probabilities of pp or nn SRCs in the nucleus are at least a factor of six smaller than that of pn 
SRCs. Our result is the first estimate of the isospin structure of NN-SRCs in nuclei, and may have 
important implication for modeling the equation of state of asymmetric nuclear matter. 
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Studies of short-range nucleon correlations (SRCs) 
in nuclei are important for understanding the short- 
distance and large-momentum properties of nuclear 
ground state wave functions. The relevant distances in 
two-nucleon (NN)-SRCs are expected to be comparable 
to that in neutron stars corresponding to 4-10 times the 
central density of nuclei [lj . Thus SRC studies are essen- 
tial in understanding the structure of cold dense nuclear 
matter. In this context the isospin content of SRCs (i.e. 
pn vs. pp and nn pairs) is important for understanding 
the structure of nuclear matter made of either protons or 
neutrons. Studies of SRCs also give the best hope of un- 
derstanding the nature of the short-range NN repulsion. 

SRCs in nuclei have been actively investigated for 
over three decades (see e.g.Q). However, experimen- 
tal studies of the microscopic structure of SRCs were 
largely restricted due to moderate momentum-transfer 
kinematics in which it is difficult to resolve SRCs. Re- 
cently, several experiments J3|, 0, la, lI Q made noticeable 
progress in understanding dynamical aspects of SRCs. 
For Q 2 > 1 GeV 2 , Refs 01a] observed Bjorken xb scal- 
ing for ratios of inclusive (e,e') cross sections of nuclei 
A to the 3 He nucleus when ig > 1.4. This confirms 
the earlier observation of scaling for nucleus-to-deuteron 
cross section ratios J3, El , and indicates directly that the 
electrons probe high-momentum bound nucleons coming 
from local sources in nuclei (i.e. SRCs) with properties 
generally independent of the non-correlated residual nu- 
cleus. 

Based on the NN-SRC picture, which is expected to 
dominate the internal momentum range of ~ 250 — 
600 McV/c, one predicts a strong directional (back- 
to-back) correlation between the struck nucleon and 
its spectator in the SRC. Experiments 0, 0, mea- 
sured triple-coincidence events for the 3 He(e, e'pp)X and 
12 C(p,ppn)X reactions, and clearly demonstrated the ex- 
istence of such directional correlations. They also re- 



vealed a noticeable momentum distribution of the center 
of mass (cm.) of the NN-SRCs. 

In this work we extend the NN-SRC model used in 
the analyses of A(p,pp)X datajlCJ . to incorporate the 
effects of the cm. motion of SRCs. This allows us to 
estimate the probability for correlated neutron emission 
following removal of a fast proton from the nucleus in 
(p,ppn) reactions. Based on this model we extract from 
the data an upper limit to the relative probabilities of pp 
and nn vs pn SRCs in 12 C. 

The measurements of 12 C(p,ppn)X reactions j(| Q 
were performed with the EVA spectrometer at the AGS 
accelerator at Brookhaven National Laboratory [ill fl2l | . 
EVA consists of a 0.8 T superconducting solenoid, 3.3 m 
long and 2 m in diameter. The 5.9 — 9.0 GeV/c pro- 
ton beam was incident along the central axis. Coinci- 
dent pairs of high transverse- momentum protons were de- 
tected with four concentric cylinders of straw tube cham- 
bers. The experimental kinematics arc discussed in more 
details later. Neutrons were detected in coincidence with 
the quasi-elastic knockout of protons from 12 C The large 
momentum transfers — t > 6 GeV 2 in these processes 
greatly improve the resolving power of the probe and 
validate the instantaneous approximation for description 
of the removal of fast bound proton in the pp — > pp sub- 
process. For each (p,pp) event, the momentum of the 
struck proton pi before the reaction was reconstructed 
and compared (event by event) with the measured coin- 
cident neutron momentum p n . Due to the ~ s~ 10 depen- 
dence of the underlying hard pp — > pp cross section, the 
scattering takes place preferentially off a bound proton 
with larg e \p2\ in the direction of the beam (minimiz- 
ing s)[Rj, and hence should lead to a significant rate of 
emission of backward correlated nucleons due to scatter- 
ing off NN-SRCs. Data confirming these characteristics 
of A(p,ppn)X reactions are shown in Fig. 1 for 12 C. The 
data show no directional correlation for neutrons with 



\p n \ below the Fermi sea level (kp = 220 MeV/c). Above 
kp a strong back-to-back directional correlation between 
f>2 and p„ is evident. 
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FIG. 1: The correlation between p„ and its direction 7 rel- 
ative to pi- Data labeled by 94 and 98 are from Ref.Q and 
respectively. The momenta are the beam momenta. The 
dotted vertical line corresponds to &f = 220 MeV/c. 

In Ref.pJ the large value of the following ratio: 



Number of (p,ppn) events (p2,Pn > kp) 
Number of (p,pp) events (p 2 > &f) 



(1) 



was extracted, which indicates that in the 250 — 
550 MeV/c region NN-SRCs must be the major source of 
nucleons in nuclei. Within the SRC model, the numer- 
ator of F is due to scattering off the pn-SRC while the 
denominator is due to scattering off any possible con- 
figuration "pX" which contains a proton with |p 2 | ~ 
250 — 550 MeV/c. All configurations such as pn- and pp- 
SRCs, the high-momentum tail of the mean-field proton 
distribution, three and more nucleon SRCs, as well as 
SRCs containing non-nucleonic degrees of freedom could 
contribute to "pX" . We define P pn /pX as the relative 
probability of finding a pn-SRC in the "pX" configura- 
tion. In this work, using a theoretical description of the 
A(p,pp)X and A(j>,ppn)X reactions, we evaluate P pn / P x 
from the above measured ratio F. 

The present theoretical description of the A(j>,ppn)X 
reaction, in which the hard pp — > pp subprocess is ac- 
companied by the emission of a recoil (kp < p n < 
550 MeV/c) neutron, is based on the light-cone distorted- 
wave impulse approximation (LC-DWIA). This is an ap- 
propriate approximation for high-momentum-transfer re- 
actions aimed at studies of the ground state properties 
of nuclei since in this case the LC momentum fraction 
a (defined below) of the nuclear constituents is approx- 
imately conserved without much distortion due to soft 
initial and final state interactions in the reaction [h3.ll5l| . 
This validates the factorization of the hard pp — > pp scat- 
tering from the soft reinteractions, allowing to express 
the A(p 1 ppn)X cross section through the product of the 
pp — > pp scattering cross section off the bound proton, 
da/dt, and the nuclear decay function, D pn , as follows: 



T p,ppn _ SP K P 
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where a p ' ppn = d 9 a/^^-^d 2 p tn and K = 
■^Vs 2 — im 2 s. Here pi, P3 and pi are incoming, scat- 
tered and knocked-out proton four-momenta and p 2 = 
P3+P4-P1- Also, s = (p 3 +p 4 ) 2 and t = (pi-p 3 ) 2 . The 
function D pn represents the joint probability of finding 
a proton in the nucleus with a = A-g^ — ^- and trans- 
verse momentum p t , and a recoil neutron in the residual 
A — 1 nucleus with a n and pt„. Pr+ = Er — p\ is the 
LC "plus" component of the A—\ residual nuclear state. 
Here, z || pi. T ppn is the nuclear transparency for the 
yield of two fast protons and recoil neutron. Eq. J3J) cor- 
rectly reproduced the angular correlation of Fig. 1 as well 
as average values of (a n ) and (005(7)) [g, UN- 

Within the LC-DWIA, the A(p,pp)X reaction is ex- 
pressed through the spectral function, S p , as follows [lCJ: 



t p,pp 



V*— (M) 2ffl(a '* ,P * +) r, 

*—' dt a 



ni' 



(3) 



where a p ' pp = d 6 a/^-^-, and T pp represents the nu- 
clear transparency for the yield of two fast protons from 
the nucleus. For SRCs we have: 

S p (a,p u P R+ ) = J2 I D ps (a,pt,a s ,pt s ,P R+ )—d 2 pt s (A) 

where the summation is over the possible types of the 
recoil particles s from the SRC. The spectral function S p 
together with the mean-field contribution to the spec- 
tral function is normalized to unity. At small inter- 
nal momenta, the LC-DWIA reduces smoothly to its 
nonrclativistic counterpart. Reasonable agreement is 
observed |10| in an extensive comparison of calculations 
based on Eq.© with measured a and pt spectra j(|. 
The ratio F defined in Eq.QJ now can be written as: 

Q max p ma* Q m<" p rr - ~ 



F 
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Pt n Ptn 



a P,Ppnd^ d 2 pt d^ d 2 ptnd p R+ 



'Pi 



(5) 



J J a P,PP^ d 2 ptd p R 



where the limits for integration are defined by the exper- 
imental conditions. The kinematics of the experiment jy, 
Q lead to the integration by Pr+ that covers all the range 
relevant to quasiclastic scattering. 

Using Eq.JSJ and relations JU,© and (0J, one can 
relate the above defined P pn / P x averaged over the mea- 
sured range of p~2, p n and P R+ to the ratio- F as follows.: 



P 



pn/pX 



F 

tJi' 



(6) 



Here T n accounts for the attenuation of the neutron and 



R = 



^I Pt f ^I "J D p ^d 2 p t ^d 2 Ptn dP m 
"""" PT tn ""'" Pte" 

"I ^J S pn ^d 2 Pt dP R+ 



(7) 



where S pn is the part of the spectral function S p related 
to pn-SRCs only. 

In the derivation of Eq.® we used T ppn « T pp T n , 
which is justified for kinematics in which two energetic, 
> 3 GcV/c protons are produced in the projectile frag- 
mentation region while the neutron is detected with 
\p n \ < 550 MeV/c in the backward direction. The accu- 
racy of this approximation is proportional to the relative 
yield of recoil neutrons due to rescattering of the ener- 
getic protons off uncorrelated neutrons in the nucleus, 
which amounts ~ 5% for the considered kinematics. 

Note that Eq.@) gives a well-defined upper limit for 
Ppn/pX ; namely, if the nucleus is transparent to the recoil 
neutron (T n = 1), all the strength of the SRC is due to 
pn-SRCs, and if the kinematic cuts cover all the domain 
relevant to NN-SRC (F = R = I) then P pn/pX = 1. 

The function R can be estimated in the NN-SRC model 
in which the decay function, D pn is a convolution of two 
density matrices representing the relative (jJojiq) and 
cm. ({^"m ) momentum distributions of the pn-SRC: 



D pn = f^ RC (Oirel,Pt,rel) ■ fi^ n m .( a c-m.,Pt,c.m.) X 

O^C. 771. 



6 [Pi 



R+ 



™ 2 + pin M\_ 2 + p\ 



ma n m{A — a c . m .) 



, (8) 



where a re i 
and pt.c.m. 



-,Pt,rel = Pt~ -^T^Ptn, Ctc.rn. = CX + 0> n , 

Pt+Ptn- Within the SRC model 0, p^ n RC 



is related to the LC density matrix of the deuteron as: 



PSB,c{a,p t ) = v(^)?^Vm 2 + fc 2 (9) 

l — a 

where ^u(k) is the deuteron wave function, and k — 

a(2-a) ~ m ' 2 (0 < a < 2). The parameter a pn (A) 
is the probability (relative to the deuteron) of having 
a pn SRC pair in nucleus A. Note that it cancels out in 
-R, while P pn / P x contains information on a pn /aNN, with 
cinn related to the probability of NN-SRCs. 

The cm. motion of the SRC relative to the (A — 2) 
spectator s yste m is described by a Gaussian ansatz sim- 
ilar to Ref.[0| with a being a parameter. This distribu- 
tion can be expressed through the LC momentum of the 
cm. of the SRC as follows: 



Pc.m.{0i,p t ) = 2m 



1 



27TCT 2 



'{2~ a )*+ P f 



(10) 



It is normalized as J p c .m.{ a TPt) < — d 2 p t = 1. 

The factorization of Eq.(JBJ) is specific to the NN-SRC 
model, in which it is assumed that the singular character 
of the NN potential at short distances defines the main 
structure of the nucleon momentum distribution, and 
that it is weakly affected by the interaction of the NN- 
SRC with the A— 2 nuclear system. This approach give s 
a good description of the spectral function of 3 i?e[l7l|. 
as well as medium nuclei |ly, with fitted values of total 
probabilities for NN-SRC which are up to 20% less than 
the values extracted in Ref. [4|, y, Q • 




FIG. 2: R as a function of a. The dashed line is for the 
same cuts applied to both the (p,pp) and (p,ppn) reactions; 
the solid line is for the cuts of the EVA experiment. 

Figure 2 shows R as a function of a, calculated within 
the above approach, as a dashed line. Here we used the 
same 0.6 < a < 1.1 and Uf < Pt < 0.55 GeV/c cuts 
in both (p,pp) and (p,ppn) reactions. It can be seen 
from Fig. 2 that R approaches unity at a — > which is 
consistent with the prediction of the SRC model with the 
cm. at rest. Figure 2 shows a large sensitivity of R to 
the width of the cm. momentum distribution. 

To extract P pn /pX from the data of Ref. [a, using 
Eq.JSJ) one needs estimates of R, T n , and F, for the kine- 
matic cuts of the experiment, which are: 

struck proton: 0.6 < a < 1.1; p 2 > p mm 

recoil neutron: 0.9 < a n < 1.4; p mm < p n < 0.55 
72° < 9 n < 132°. (11) 

Here we follow Ref.@ and set p mm = 0.275 GeV/c to 
suppress mean field contribution to < 1% 

R: The solid curve in Fig. 2 represents R calculated 
for the experimental cuts of Ea. l|ll|l . It is different from 
unity even for a — > due to the different integration 
ranges for p2 and p n ; Ea. Hll|) places no restriction on 
p 2t = Pt , while Pnt > 204 MeV/c. 

T n : We used T n = 0.85 obtained from the simulation 
based on Eq. J2J) using the experimental limits of Eq. pill . 
T n includes Pauli Blocking effects according to Ref.[la| 
as well as a small correction due to neutrons from the 
rescattering of the fast incoming and outg oing protons 
calculated in the eikonal approximation [l9t |20| | . 

F: The ratio F is taken from Ref.jTfl , where it is quoted 
as F = 0.49 ± 0.13 for p min = 0.22 GeV/c. For p min = 
0.275 GeV/c our analysis of the same data gives: 

F = 0.43±8#- (12) 

The uncertainty is dominated by the low statistics. Since 



the lower limit of the extracted P, 



pn/pX 



is sensitive to the 



lower limit on F (Eq.©) we choose for F the best one 
standard deviation low limit allowed by the data. 

Using the above values of R, T„, and F, we estimate 
Ppn/pX from Eq.©. Figure 3 shows the a dependence 
of Pp n /px for F = 0.36, 0.43 and 0.55, respectively. 
Since P pn / P x < E there is an interesting correlation 
between a and P pn /pX, which allows us to put a con- 
straint on a . For example for F = 0.43, a can not exceed 
174 MeV/c. To evaluate P pn / P x we use the magnitude 
of a exp = 143 ± 17 MeV/c extracted from the same data 
set 0- This value is in excellent agreement with the 
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FIG. 3: Pp n /px as a function of a. The shaded area corre- 
sponds to the Pp n / P x values at a exp — 143 ± 17 MeV/c. 



theoretical expectation of 139 MeV/c of Ref. 16] . Note 
that a exp dictates the range of possible values for P pn / P x ■ 
From the central (0.43) and minimal (0.36) values of F 
we obtain: 



P, 



pn/pX 



Q2+ 008 

u - yz -0.18- 



(13) 



This result indicates that at least 74% of the time the 
removal of a fast proton is accompanied by the emission 
of a fast recoil neutron. It allows us also to estimate an 
upper limit of the ratio of absolute probabilities of pp— 
to pn— SRCs (p^), by assuming that pp-SRCs carry all 
the remaining strength in "pX" configuration|2l|: 



P, 



04+ 009 



(14) 



P 1 

-ET^ ^ o( 1_ fpn/ P X> 
r pn ^ 

This result can be used to estimate separately the ab- 
solute probabilities of pn, pp and nn SRCs in the nu- 
clear wave function. For this we use the total proba- 
bility of NN-SRCs (Pnn( 12 C) = 0.20 ± 0.042) obtained 



2 and xb inclusive 
The value of Pnn 



by combining the results of large 

A(e,e')X data from Refs.0, EM 

is practically independent of whether or not we account 

for the cm. motion of the SRC. Using Pnn together 

with Eq.lO we obtain P pn ( 12 C) = 0.184 ± 0.045 and 

P pp ( 12 C)=P nn ( 12 C)<0m. 



Summarizing, from Eas. H13ll4|) we conclude that our 
analysis of (p, pp) and (p, ppn) reactions indicates that if 
a nucleoli with momentum between 275 — 550 MeV/c is 
removed from the nucleus using a high momentum and 
energy transfer probe, at least 74% of the time it will 
originate from pn-SRCs. The data also show significantly 
smaller (a factor of 6 at least) probabilities for pp and nn 
NN-SRCs as compared to pn NN-SRCs in Carbon. This 
result may have an important implications for modeling 
the equation of state of asymmetric nuclear matter. It 
indicates that the equation of state of neutron stars could 
be affected by the large suppression of the number of neu- 
trons above the Fermi sea. Future experiments involving 
high-Q 2 (e, e'pp) and (e, e'pri) reactions off medium-to- 
heavy nuclei (e.g., Ref.J22|), as well as theoretical studies 
based upon realistic treatments of SRCs and final-state- 
interaction effects (see, e.g., Ref.|23ll24l l25|L will allow 
one to verify and significantly improve the accuracy of 
the present results. 
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